HIGHER ORDER DERIVATIVE ESTIMATES FOR 
FINITE-DIFFERENCE SCHEMES 



ISTVAN GYONGY AND NICOLAI KRYLOV 

Abstract. Wc give sufficient conditions under which solutions of 
finite-difference schemes in the space variable for second order pos- 
sibly degenerate parabolic and elliptic equations admit estimates of 
spatial derivatives up to any given order independent of the mesh 
size. 



1. Introduction 

This is the second part of a series of papers devoted to studying 
the smoothness of solutions to finite difference schemes for parabohc 
and elhptic partial differential equations given on the whole M"'. These 
equations can degenerate, for example be just first order PDEs. As in 
[Hj, the first part of this series, we consider a grid in M.'^ and a large 
class of monotone finite difference schemes in the space variable x in 

For each small parameter h > the given grid is dilated by h and 
for each a; G M'' it is shifted so that x becomes a mesh point. We are 
interested in the smoothness in x of the solution Uh{t,x) of the dif- 
ference scheme. In [13] estimates, independent of h, for the first order 
derivatives of Uh in x were obtained under general conditions introduced 
there. In the present paper we investigate the higher order derivatives 
of Uh in X. The main results give estimates, independent of h, for the 
derivatives of Uh in x up to any given order m. The conditions extend 
those from [H]. Using these results in the continuation of this paper 
we estimate the derivatives of Uh in h, and that allows us to develop a 
new method of obtaining the power series of Uh in h. Hence we get ac- 
celerated finite-difference schemes by using Richardson's extrapolation. 
Namely, under general conditions we show that the accuracy of finite 
difference schemes for parabolic and elliptic PDEs can be improved to 
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any order by taking suitable linear combinations of finite difference ap- 
proximations with different mesh-sizes. For elliptic PDEs this result is 
announced by Theorem 2.4 in [T3] and for parabolic PDEs by Theorem 
12.31 below. We hope to develop these results in domains for uniformly 
nondegenerate equations later. 

Derivative estimates for finite-difference approximations for linear 
and for nonlinear PDEs play the paramount role in establishing the 
rate of convergence of the approximations. The importance of such 
estimates is demonstrated recently by [2U], [2T] and [22], presenting 
the first rate of convergence result in the sup norm of finite-difference 
approximations for fully nonlinear degenerate Bellman equations. Ideas 
from these publications are used and developed further in [2], [17], [3], 
[Z], [H], [9] and [16]. Recent results on estimating the Lipschitz constant 
and second order differences of finite-difference approximations for a 
large class of fully nonlinear degenerate PDEs, including the normalized 
Bellman equations are presented in [21]. In [8] first order derivatives 
of finite-difference approximations to degenerate linear parabolic and 
elliptic PDEs are estimated and are used to establish sharp estimates 
on the rate of convergence of the approximations in the sup norm. 

Finite-difference methods for solving PDEs have been extensively 
studied since the first half of the last century. Let us mention the pio- 
neering papers by R. Courant, K.O. Friedrichs and H. Lewy [6], S. Ger- 
schgorin [12], and publications by D.G. Aronson, J. Douglas, F. John, 
H.O. Kreiss, O. Ladyzhenskaya, P.D. Lax, W. Littman, Lyusternik, 
J. von Neumann, I.G. Petrovskii, A. A. Samarskii, G. Strang, A.N. 
Tikhonov, V. Thomee, O.B. Widlund and many others (see, e.g., [T], 
HO], [IB], n 23, EH], [21], [30], [32], [33], [10], jH] and the refer- 
ences there.) We refer also to the review paper [3S], handbook [30], 
and well-known monographs and textbooks for more information on 
the subject ([5], [11], [13], [20], [31], [31], [35], [33). 

The paper is organized as follows. The main results. Theorems 12.11 
and 12.21 are presented in Section [21 Here we formulate also a result. 
Theorem 12.31 on accelerated finite difference schemes, which we will 
prove in the continuation of this paper by using Theorem 12. 1[ As we 
have pointed out above, the idea of the proof of Theorem 12.31 is based 
on a power expansion of in h. This idea was already applied by 
the authors to show how to accelerate other approximation schemes 
(see, for instance, [15]) and it seems to the authors that it was never 
used before in the framework of finite difference schemes in the sup 
norm for degenerate elliptic and parabolic equations although much 
effort was applied to developing this and other methods of improved 
approximation for uniformly nondegenerate equations in domains (see. 
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for instance [1] and the references therein). It is worth saying that, in 
contrast with [1] and many other papers deahng with the expansion, 
we do not use any information from the theory of PDE and, as a 
matter of fact, the existence of smooth solutions for degenerate eUiptic 
and parabohc equations follows directly from our results. We deduce 
Theorem 12.21 from Theorem 12. and conclude Section [2] with verifying 
the rather delicate conditions of the main results. Assumptions 12.41 and 
12.51 for a class of examples, given in Remark 12.21 before the formulation 
of the theorems. The proof of Theorem 12. II is given in Section[3l and the 
final section. Section HJ is devoted to further discussions of Assumptions 
[231 and [231 



2. Formulation of the main results 

We take some numbers ho, T G (0, oo) and for each number h G 
(0, /lo] "we consider the integral equation 

u{t,x) = Qhix) + {Lhu{s,x) + fhis,x)) ds, {t,x)eHT (2.1) 
Jo 

for u, where g = gh = gh{x) and f = fh = fh{s, x) are given real-valued 
Borel functions of x G M"^ and (s, x) G Hx = [0,T] x M*^, respectively, 
and L = is a linear operator defined by 

Lh^{t, x) = L^Lfit, x) - c{t, x)lp{x), 

L^^{t,x) =Ll^{t,x) = j^Yl (lxit,x)^h,\^ix) + '^pxit,x)6h,x^{x), 

AeAi AeAi 

for functions (f on M*^. Here Ai is a finite subset of M."^ such that 
^ Ai, and px{t,x), qx(t,x) are real-valued functions of {t,x) G i^oo = 
[0, oo) X M*^ given for each A G Ai, and 

6x(p{x) = Sh,xf{x) = l-{(p{x + h\) - fix)), A G Ai. 

n 

As usual, for multi-indices a = {ai, . . . ad), ai = 0,1,..., we use the 
notation 

* i 

For smooth ip and integers A; > we introduce D'^ip as the collection 
of partial derivatives of (p of order k, and define 



4 



I. GYONGY AND N. KRYLOV 



Let m > be a fixed integer and let Ki G [1, oo) be a constant. We 
make the following assumptions. 

Assumption 2.1. For any A G Ai the derivatives in x of px,qx, c, /, g 
up to order m are continuous functions in (t, x) G Ht and, for k = 
0, ...,m and some constants we have 

sup ( J2 (\D'qx\' + \D'px\') + \D''c\') < Ml (2.2) 

By Theorem 2.3 of [2] under Assumption 12.11 for each h G (0, /iq], 
there exists a unique bounded solution Uh of (12.1 p and this solution 
is continuous in along with all its derivatives in x up to order m. 
However, the bounds, provided by this theorem for these derivatives 
depend on the parameter h. Our aim is to show the existence of bounds, 
independent of h, if in addition to Assumption 12.11 the assumptions 
below also hold. 

Assumption 2.2. For alH G [0, T] 

'^^\qx{t,x) is independent of x. (2.3) 

AgAi 

Introduce 

X\ = Xh,x = q\ + hpx- 

Assumption 2.3. For all {t,x) G Ht, h G (0, Hq], and A G Ai, 

Xx{t,x)>0. (2.4) 

There exists a constant cq > such that c> cq. 

Obviously Assumption 12.31 implies that qx > 0. 

Remark 2.1. The above assumption: c > Cq > 0, is almost irrelevant 
if we only consider (12.11) on a finite time interval. Indeed, if c is just 
bounded, say \c\ < C = const, by introducing a new function v{t,x) = 
u{t, x)e~^'"* we will have an equation for v similar to (12. ip with L^v — 
(c + 2C)v and /e"^*"* in place of Lu and /, respectively. Now for the 
new c we have c + 2C > C. 

Take a function tx defined on Ai taking values in [0, oo) and for 
A G Ai introduce the operators 

Tav? = Th,xLp{x) = ip{x + hX), 6x = 5h,x = Txh~^{Tx - 1). 

For uniformity of notation we also introduce A2 as the set of fixed 
distinct vectors none of which is in Ai and define 

= ^h,t = To-Dj, Tgi = Tfi^ii = 1, A = Ai U A2, 
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where tq is a fixed parameter satisfying 

To > 0. 

For integers k = 1,2, ... and A* G A, i = 1, 2, k, introduce tlie multi- 
vectors 

A = (A\...,A^') G A'^ 

and tlie operators 

Tx = Thx = Thxi...Tf^y^k, 5x = Sh,x = Sh,x^---Sh,x'^- 

It is also convenient to set A^ = Ag = A° = {0} and to introduce 
6o = 6q and To as unit operators. For /i G A'^ and k < m we set 

Qip = ^ qxSx^, Qt,^ = {S^qx)5x^, 

AeAi AeAi 

^ Pa^aV^, P,,^ = ^ {6^px)5x^, 

AeAi AgAi 

A,{cp) = 2 J2 i^x^)LlTxip, Q(^) = Yl X.iS.iff. 
agA'= MeAi 
Below B{M.'^) is the set of bounded Borel functions on M"^ and ^ is the 
set of bounded operators K, = K,h = }Ch{t) mapping B{M.'^) into itself 
preserving the cone of nonnegative functions and satisfying ATI < 1. 
Set 

\A^\'=J2\M^ \\Mr=J2\rxX\^ 

AeAi AeAi 

Finally, fix a constant 6 G (0, 1]. 

Assumption 2.4. We have m > 1 and for any h G (0, ho], there exists 
an operator /C = /Ch,m G .ft, such that 

mAii^) < {l-6)J2Qi^xv)+KMv) + 2{l-6)c}C{J2\hfn (2.5) 

AeA AeA 

on Ht for all smooth functions ip. 

Assumption 2.5. We have m > 2 and, for any h G (0, /iq] and n = 
1, ...,171, there exists an operator /C = JCh^n ^ -ft, such that 

nJ2MM+n{n-l) J] (M^aTa^^ < (1 - 5) Q(5a¥^) 
i^eA AeA2 AeA2 

+ KiJ2 2(^A<^) + 2(1 - 5)c/C( ^ |5a<^|2) +K^IC{J2 (2-6) 
AeA AeA2 AeA 

on Ht for all smooth functions (f. 
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Obviously Assumptions 12.41 and 12.51 are satisfied if q\ and p\ are 
independent of x. In the general case, as it is discussed in [13], the 
above assumptions impose not only analytical conditions, but they are 
related also to some structural conditions, which can somewhat easier 
be analyzed under the following symmetry condition: 

(S) Ai = — Ai and q\ = g_A for all A G Ai. 

Notice that, if condition (S) holds then 

h-^ J2 qxit,x)6x^ix) = (1/2) J2 qxit,x)AMx), 
AeAi AeAi 

where = A/^ a and 

^{x + h\) - 2<^{x) + y?(x - h\) . . ( ^ 
Ah^x^{x) = — = -d-xOx^{x). 

Remark 2.2. In this remark we suppose that Assumptions 12. ![ 12. 2[ 
and 12.31 hold and m > 2. Assumption 12.41 is discussed at length and 
in many details in [14] . At the end of this section we show that if 
condition (S) holds and , for all A G Ai, tx > and qx > k, where 
K > is a constant, then both Assumptions 12.41 and 12.51 are satisfied 
for any Cq > and S G (0, 1), if ho is sufficiently small and Tq > 0, Ki, 
and /C are chosen appropriately. It follows from Remark 6.4 of [H] and 
Remark 14.21 that the above condition k> can be dropped, provided, 
additionally, that cq is large enough (this time we need not assume that 
h is small). 

By the way, as we have seen in Remark |2. 11 the condition that cq be 
large is, actually, harmless as long as we are concerned with equations 
on a finite time interval. 

Mixed situations, when c is large at those points where some of qx 
can vanish are considered in Section m 

Now we are in the position to formulate our main results. 

Theorem 2.1. Let Assumptions \2.1\ through \KR hold. Then for h G 
(0, ho] we have 

m 

sup l^'^^l < N^Fm + G^), (2.7) 



where 



Fn = Ysup\D''fh\, Gn = Y sup \D''gh\ 



k<n k<n 

and N depends only on To, m, 5, Cq, Ki, |Ai|, ||Ai||, Mo,...,Mm. 
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We prove this theorem in Section [3l Now we derive from it an esti- 
mate for the solution of the equation 

LhV + fh = inM^ (2.8) 
when qx, px, c, and / are independent of t. 

Theorem 2.2. Let Assumptions \2. 1\ through \2.4\ be satisfied. Suppose 
that qx, Px, c, and f are independent of t. Then the following state- 
ments hold. 

(i) There exists a unique bounded solution v = Vh[x) of (12.81) . More- 
over, all derivatives in x of v up to order m are bounded continuous 
functions on M.'^. 

(a) Let Assumption \2.5\ be also satisfied. Then 

m 

sup < ATF™, (2.9) 

where 



Fm = ^supl^V/il 



fc=0 

and N depends only onr^, m, 5, Cq, Ki, \Ai\, \\Ai\\, Mq, M^. 

Proof. Statement (i) is proved in [T^ (see Theorem 2.2 there). To 
prove (ii) take u = CqJ, where 7 > is so small that c — v > cq/2 and 
conditions (12. 5p and (12. 6p hold with c — v and 5/2 in place of c and 5, 
respectively. Define u{t,x) := v{x)e'^^ and observe that u satisfies 

= L\ - (c - u)u + e^'f. 
dt \ J J 



By Theorem EH] for x G M"* 

m m 

e''^^|D'^w(a;)| = ^ |Z^'=n(T, x)| < A^e^^F^ + ^ sup |D\;(a;)|. 

k=0 k=0 k<m 

By multiplying the extreme terms by e~^^ and letting T ^ 00, we get 
the result. □ 

The above theorems have important applications in the numerical 
analysis of finite difference schemes for parabolic and elliptic PDEs. 
Using them in the continuation of the present paper we obtain accel- 
erated finite difference schemes for second order (possibly) degenerate 
parabolic and also for second order (possibly) degenerate elliptic PDEs. 
In particular, we will consider the Cauchy problem 
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d 

—u{t,x)=Cu{t,x) + fo{t,x), te(0,T],xeM^ (2.10) 

u{0,x) = go{x), xeR"^ (2.11) 

with the operator 

d d 

£ := i ^ ^ qxKXjDiDj + ^ ^paAjA - c. 

AgAi i,j=l AgAi i=l 

By a solution of fl2.10p - fl2.lip we mean a continuous function u(t, x) on 
i^T, such that for each t it is twice continuously differentiable in x, is 
bounded in Ht along with its derivatives in x up to second order and 
satisfies ^ 

u{t,x) = go{x) + / [Cu{s,x) + fo{s,x)]ds 



in Ht- 

To formulate one of the main results of the continuation of the paper 
we fix an integer k >0 and set fh = fo,gh = go, 

k 

= ^ bjU2-jh, 
3=0 

where M2-Jh is the solution to (12. ip with 2~^h in place of /i, 

{h,M.-M) := (i,o,o,...,o)i^-\ 
and is the inverse of the Vandermonde matrix with entries 
yij 2-('-i)0--i), i^j = 1,...,A; + 1. 

Theorem 2.3. Let Assumptions \2.'J\ and \2.1\ with m > 3{k + 1) hold. 
Also let condition (S) be satisfied. Then fl2.10l) - fl2.1ip has a unique 
solution uq, and 

\uh{t, x) - uo{t, x) I < Nh^+^ (2.12) 

holds for all [t, x) G Ht, h G (0, ho], where N is a constant depending 
only on T , k, d, |Ai|, Hq, the number of elements in Ai, on Mq, . . . .Mm, 
on sup^gfo^r] l/o(i)|m and on \gQ\m- 

Now we prove Remark 12.21 Instead of condition (S) we assume the 
following weaker condition 

(S'): Ai = -Ai and Dqx = Dq_x for A G Ai, 

and proceed with the proof of Remark 12.21 as follows. Clearly, 



5^(5a^)L°TaV9 = /i + /2, 



aga 
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with 

h ■■= (^A^)L°T,^, h := (^A^)I^^^. 
AeAi AeA2 
Due to condition (S') 

AgAi AgAi 
A,fiGAi A,/^gAi 

+ Y CSxv)CSxX,)5,5xV =■■ ^ + I?^ + 

A./iGAi 

T _ r(l) , r(2) 

I2 — I2 + ^2 ^ 

where in the notation C, = DLp/\D{p\ and = C,iDitp, 

d 

= ^0 Y E iD,^mq,)A,v = |ro^|/^^| Y %ii)^>^^^ 

j=l aiGAi ^GAi 
d 

j=l AiGAi /xGAi 

Set 

r := maxTA, r := min Tx, 

AgAi AgAi 

and observe that xx > i^f^ > for sufficiently small h, and that 
c > Co > 0. Then by Young's inequality, we obtain 

AgAi 

AiGAi AGA2 

2m/f <roiVcolc^|5A^|^ 
aga 

where is a constant depending only on m, k, 6, r, r, the number of 
elements in Ai, and on the supremum norm of the gradients of px and 
qx in X. Summing up these inequalities and taking tq > sufficiently 
small we get ( I2.5P with Ki = 3N, unit operator /C, and with 6 as close 
to 1 as we wish. 

This result is obviously applicable to 6,yip is place of ip for any i/ G A. 
It follows that for any S G (0, 1) and appropriate constant Ki we have 

mYMM<a-S)YQCSxV^ 

i^GA AgA2 
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+ KrY, 2(^AV') + 2(1 - 5)c V^M^- (2-13) 

AgA AgA2 

Now we show that Assumption 12.51 holds. Clearly, 

^ {5xv)Q\Tx^ = ^ {5x'^)QxTx'^ + 2 ^ {5x'f)QxTx'f 

AeA2 AgA2 AGAixAa 

+ J2 (^AV')Qa^ = Ji + /2 + h. 
AGAi 

Using 

A^Ta = A^. + h5xi6x2i5^ + 6.^) + (5^1 + ^aOI'^m + ^-m) 
for A = (A\ A^) G and yU G Ai, we have 

AGAf 

with 

A^'^ = I E {h^)Chq,){^Sx^-S.,)6,ip, 
aga2,^gAi 

As above, we have 

AGAf A,/xgAi AgA 

E 

aga2,^gAi agA2 
and hence, by Young's inequality 

AgA 

n{n - ) < h^N Y Qi^^f) + A^E (2.15) 

AGA2 AgA 

where A^ is a constant depending only on m, n, f, the number of 
elements in Ai and on the supremum norm of l^^gAl- Similarly, 

d 

n{n - l)/2 = n{n - 1)tq ^ ^ (6uToDi(p){6uDiqf,)T^A^ip 

i=l 1/,/iGAi 

< roA^E + ^^oclC{J2 ^z^*^)' (2.16) 

AgA /xgA2 
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d 

n{n - l)h = ln{n - 1)t^ ^ ^ {T^Dij<^){Dijqf,)Af,<^ 

i,j=l fJ.eAi 

< Nrlc J2 (2.17) 

AGA2 

where N denote some constants depending on m, d, Cq, k, f, r, the 
number of elements in Ai and on the supremum norm of \D'^q\\. Sum- 
ming up the inequahties fl2.13p through (12.171) and choosing tq and ho 
sufficiently small we obtain fl2.6p . 

3. Proof of Theorem 12.11 

For m = 1 estimate fl2.9l) holds by virtue of Theorem 2.1 from jTi] . 
proved by the aid of the following version of the maximum principle 
(Corollary 3.2 in [Mj). 

Lemma 3.1. Let Assu'mption \2. 1\ with m = be satisfied and let Xa ^ 
for all A G Ai. Let v be a bounded function on Ht, such that the partial 
derivative DtV := dv(t,x)/dt exists in H^- Let F be a nonnegative in- 
tegrable function on [0,T], and let C be a nonnegative bounded function 
on Ht such that 

V := sup(C — c) < 0. 

Ht 

Assume that for all (t, x) G Ht we have 

DtV < Lv + Cv+ + F, (3.1) 
where v{t) = snp{v(t,x) : x G W^}. Then in [0,T] we have 

v{t) < v+{0) + sup F, (3.2) 

m 

where a+ := (|a| + a)/2 for real numbers a. 

For the proof of this lemma we refer to [H]. In order to obtain 
Theorem 12.11 for m > 2 we need some more lemmas. First we prove 
a lemma which will be used a few times in the future. By /C in the 
lemma and later in the article we mean a generic operator of class ^. 
This operator may change each time it is mentioned even in one line 
(cf. the use of o{n)). Thus, for example, for nonnegative functions a, 
P on M'^ the formula a/C + PJC = (a + /3)/C means the simple fact that 
for any /Ci , /C2 £ 

a/Ci + /3/C2 = (a + /?)/C3 

with 

^3 := ^/Ci + ^/C2 G (§:=0). 
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Lemma 3.2. Let Assu'mption \2.1\ be satisfied. Let n > 1 be an integer. 

Set 

A' = \\A^r + rl 
Then for any (f G we have 



Furthermore, if 1 < n < m, then for any ip E C" 

AgA" -'^^ AtGAi 



(3.3) 



<r,'\A,\'A'-{snp I^>mI')^( E I^^V^I') ' (3.4) 

MGAi AgA 

and if assumption fl2.3p holds, 

AgA" ■'^^ AiGAi 



< -0^^' 



^|Ai|M2"(sup 5^ |/^"g,P)/C( 5^ \6M')- (3.5) 

MGAi AGA2 

Finally, 

( E ^aI^aV')^)' < J](5,v^)^ < Mo^( J](5,^)^)^ 

AgAi AgAi AgAi 

Proof. It is easy to see that for A G A" we have 



where, for y{X, 9) = \^6^ + ... + \^9^ and Tx = txiTx^ ■ ... ■ r^n, 

d 

M9,x)=Tx E ^i---K„D,,-...-D,M^ + y{\,9)). 

il,...,in = l 

By Cauchy's inequahty 



\6x^ix)\' < / \M9,x)\'d9, 

J[0,h]" 

\M0.x)\<Tx\\\\DM^ + y{^,0))\, 

where |A| := |A^| ■ ... ■ jA"']. It follows that the first inequality in (13.31) 
holds with with Vh in place of /C, where 

P,^(x) = ||Ai|r2«^^2|^|2^-„, f ^^^ + y^x,9))d9 (i:=0). 

AGAf -^iOM" 
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Since 

E-^2|\|2 II A ||2n 
^x\M = ll^lll ) 

Vhl < 1, that is "Dfe e ^ and the first inequahty in fl3.3p is proved. 

To prove the second one introduce T>h^k as the operators for which 
first inequahty in fl3.3p holds with k in place of n, recall that = A2 = 
{0} and 60 is the identity operator, and observe that the left-hand side 
of the second inequality equals 

n n 

n 
k=0 

with Sh G This proves the first assertion of the lemma. 

To prove (13. 4p notice that by Cauchy's inequality and (13. 3p for A G A" 

MGAi 

< E (^^^p/^)'^^ E ('^^^)' ^ iAip E(^^^'^)'^^^^.i(i^'^i')- 

Hence the left-hand side of (13.40 is less than 

|AiP E CSxP,rTxV,,^i\D^\') 
mgAi,agA" 

=: |AiPa2"(supEI^>mI')^/.(I^^I')- 
Here Tth G ^, since by (13. 3p 

E (^apm)' < E i^vn < ^''^sup E 

This proves (13. 4p . To prove (13.50 . notice that Xl^/^^AQ'^t = 0, which 
implies that 

QxTx^ = E i^xQ^i)Txi5^^ - ^iiDiip) = E (^AgM)^A^^M' 

AtSAi mSAi 

where 

Hence as above the left-hand side of (13. 5p is less than 

E J2C5xQ.rTx E i^/^i' =■ ^"(^^p E id^'^.^M E i^/^i')' 

AeA" /^6Ai /xsAi ■'^^ /^6Ai /iSAi 
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where J^h ^ ^- Furthermore, 



ipfj,{x) = h I (h- 9)iJ,iHjDijLp{x + n9) dO 
Jo 

< / \D^^{x + fi9)\d9=\fi\'^ [ \D^^{x + hne)\d9, 



\i^t.{x)r < l/^r / \D'ip{x + M) r de, 

Jo 

and we obtain (13. 5p with the operator 
which is in ^ because 

/CI = |Ai|-^ J2 l/^l^ ^ 1- 

Since the last assertion of the lemma is obvious, the lemma is proved. 

□ 

The following lemma can be proved easily by induction on n. (Sums 
over empty sets of indices are defined to be in the lemma, and every- 
where in the article.) 

Lemma 3.3. Let n > 1 be an integer, ifj and (p be n times continuously 
differentiable functions on M*^, and A G A"-. Then 

n 
i=l 

l<i<j<n 

+ i^Kii,-,ih)'^)h{h,...,ik)T\{h,...,ik)^ 
l<i\<. . .<if^<n 

+ ... + {5xilj)Txip. (3.6) 

where X{ii, ...,ik) = (A*^, A**^), X{ii,...,ik) is the sequence of vectors 
A^,...,A" from which the vectors standing on the places with numbers 
ii, ...,ik are removed, and Sx^-^ := 1 for n = 1. 

Proof of Theorem [2lll Recall that A° = {0} and 6o = Tq is the 
unit operator. Fix h G (0, ho], for < k < m set 

u = Uh, f = fh, Vo = u'^, Vfc = |(5auP, Vk{t) = supVk{t,x), 
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and recall that F„ is introduced in Theorem I2.1[ Take an integer n G 
[1 , m] . Then we have 

LlVn = 2 iSxu)LlSxu + Qihu). (3.7) 

AeA" AeA" 



By Lemma [3.31 

AeA" AgA" n>fc>l 

where 

By Assumption | 



AgA 

+ 2(l-5)c/C(5^|5A5M«r). (3.9) 

AeA 

Hence, 

< (1 - 5) 5^ Q(M + /Ci 5^ Q(M + 2(l-5)cK. 

AeA" AeA"-i 

Next, if n > 2, then 



fteA"-i ^teA"-2 

= XI ('^5Za(4(M) + I^(^-i)^2(M)' 

/ieA"-2 !/ga 
so that by Assumption 12.51 

^jgA"-2 AgA2 

+(1 - 5) 5^ Q(5An) + i^i 5^ QiSxu) + 2(1 - 5)cK + 

AeA" AeA"-i 
By Lemma [3.21 

n-l 

n(n-l) 5^ 5^(5aW^aTa5^w+ ^ ^n,k < S'Vn + N^Vk, 

/ieA"-2 AgA2 n>fc>3 fc=l 
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where and below by the sum over an empty set we mean zero. It follows 
that, for n G [1, m], 

LlVn > 2 J2 (^An)5Alv°,n + 6 J2 2(^An) - 111 Yl 

AeA" agA" AeA"-i 

n-l 

- {2c-25c + 6'^)Vn- Nj^Vk- (3.10) 

fc=i 

Next, 

2 ^ {6xu)6xLlu = DtVn + Ri - R2, (3.11) 

AeA" 

where 

Ri:=2j2 i^xu)5xicu), R^ := 2 i^xu)Sxf. 

AeA" agA" 

Similarly to 



/2i = 2cK + 2^C„% 



A;=l 

where 

-Rlfc := E ('^M^A^)(^A'C)T^(^A^)- 

AgA"-'' /xeA'= 

By our assumptions, (13.31) . and Cauchy's inequality 

\WT,\hu\<{ Y iVhurfXi Y Chuff' 

xeA"-'' AeA"-*^ AeA"-'^ 

AgA"-'' 

|2^1/2^ ^fc^Sup I nfc^lU/l/2T7V2 
MgAfc 

Cauchy's inequality also allows us to estimate R2 and conclude from 
IKWf that, for n G [l,m], 

- 2cK - DtVn > 6Qn 
n-1 

- KiQn-i - (2c - 25c + 26^)Vn -Nj^Vk- NF^ (3.12) 

fc=0 

where 

Qk=Yl QCxu). 

xeA'' 

We now prove (12.91) by showing that for each n G [l,m\ 

Vk<N{Fl + Gl), fc = 0,l,...,n. (3.13) 
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We prove this by induction on n. By Lemma [3.11 we have 

V,<N{F',+Gl). 
Using this, from fl3.12p we obtain 

Vi < N{F^ + Gl) 

by Lemma 13.11 provided that < 2c — 25c + 26^ < 2c — S"^ which is 
true indeed if 

36 < 2co. (3.14) 

This may look hke a nontrivial restriction on 6. However, obviously, if 
our assumptions are satisfied with a. S e (0, 1), they are also satisfied 
with any 6' G [0,6]. Therefore, without losing generality we suppose 
that (13.141) is valid. Thus we have obtained (13.131) for n = 1. Let n > 2 
and assume that (13.131) holds with n — 1 in place of n. Then from (13.121) 
for k = 1, ...,n we get 

Lyk-2cVk-DtVk > 6Qk-KiQk-i-CsIk=nVn-N{F^ + Gl), (3.15) 

with Gs = 2(c — (5c + 5^). Actually, (I3.15P is true also for = if we 
set Q„i = 0, since 

Ll{u^) - 2cu^ - Dt{u'^) = 2u{Llu - cu - Dtu) + Qq{u) 

= -2uf + Qoiu) > -N{F^ + Gl) + Qoiu). 

Next we set fi = Ki/6, multiply (I3.15P by /i""^' and sum up the resulting 
inequalities with respect to A; = 0, ...,n. Then, for 

n 

we obtain 

LlWn - 2cWn - DtWn > SQn - G - N {F^ + Gl) 

> -GsWr.-N{F^ + Gl). 

Recalling (I3.14p and using Lemma 13.11 shows that (13.130 holds. This 
justifies the induction and proves the theorem. 

4. Discussion of Assumptions 12.41 and 12.51 

In [13] there are many sufficient conditions for Assumption 12.41 to be 
satisfied. In this section we suppose that only Assumptions 12. ll and 12.31 
are satisfied and m > 2. Assume also that for a number f > we have 
that, for any A G Ai, 

either T\>f or Dqx(t,x) = for all it,x). (4.1) 
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Recall that by /C we denote a generic operator from class which may 
depend on h and t, and may change each time it is mentioned even in 
one line. 

Remark 4.1. Assume that m > 2, Ai = — Ai, qx = and p\ > 0. 
Then, since qx are twice continuously differentiable in x and nonnega- 
tive by Assumption [231 '^6 know that := ^/q\ is Lipschitz continuous 
in X with the Lipschitz constant independent of t. 

Conditions (12. 5p and (12.61) involve a mixture of finite differences and 
derivatives. Therefore, it is reasonable to try to find conditions in terms 
only of finite differences which would imply (12.51) and (12. 6p . We claim 
that (12. 5p and (12. 6p are satisfied with a tq > and, perhaps, different 
6, /C, Ki if for all smooth ip on Ht and n = 1, m we have 

2m i^xV)LlTxV <{1-S)J2 2(^a¥^) 

AgAi AgAi 

+ K^Q{^) + {l-6)clC{J2f^M')^ (4-2) 

AgAi 

2n ^ {SxSu^)LITxSu^ + n{n - 1) ^ {5x^)Q\Tx^ 

A,i/GAi AgA2 

< (1 - 5) ^ QiS.if) + K,J2 
+ (1 -5)c/C( J2 Ih^n+KM J2 (4-3) 

AGAf AgAi 

(Notice that the term 2(1 — 5) in (12.51) and (12.61) is replaced now with 
1-S.) 

To prove this, first observe that as it follows from Remarks 5.1 and 
5.2 of [H], owing to (14. 2p . the above mentioned properties of r^, and 
(14.11) . there exist constants 6, tq G (0, 1], and Ki, perhaps different from 
the above ones, such that 

mA,{^) < {l-5)J2Qi^x^) + KM^) + {l-5)clC{J2\^x^n (4-4) 

aga aga 

on Ht for all smooth functions ip provided that tq G (0,fo]. Thus, a 
condition even somewhat stronger than (12. 5p is satisfied. Next, observe 
that the left-hand side of (12.61) equals 

B + A' + A'^ + a; + B' + b\ 

where B is the left-hand side of (14.31) . 

A' = n^ Ai(4</^), Al = 2n {6jx'p)QuSx¥>, 

ueA2 AGAi,i/GA2 
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A; = 2n CUxV)PJx^, 

AGAi,i^eA2 



2 



B' = 2n{n-1) ^ {6xv)QxT^^, 

AeAixA2 

B" = n{n-l)J2i^xV)QxV■ 
Here by (jOD 

A'<{l-6) Yl Qi^xSuV) + KiJ2Qi^'^y') 

ueA2,xeA ueA2 

+{l-6)clC{J2fSxv\")- 

AeA2 

Then, 

d 

A,/xeAi j=i 
d 

A,mGAi j=l ^gAi.AgAf 

AeAi,i/eA2 ^gAi,AeAf 

< (l/16)5c 5^ (6,^)' + iVro^ 

AeA2 agA2 
where and below by we denote various generic constants independent 
of ip, (t,x), and Tq. Next, quite similarly 

d 

A; = 2nT^ Y $^P,5a^)P,Pm)Ma^ 

A.MGAi j=l 

< {l/16)6cr^ Y T.^Dr5,^f + Y (^^A^)' 

AgAi 3=1 AgA2 

< (l/16)5c Y (^AV')' + NtI Y (^a^)'- 

AeA2 agA2 

Now we estimate B' and B" . We have 

d 

B' = n{n-l)rl Y $^PAv')PA?M)A^T;,y.. 
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Here 

and it is seen that 

d 

B' < (l/16)5cro^ H^DfS^^f + Nrl ^ [hv? 

AeAi i=i AeA2 

< (l/16)5c {hvf + Nrl Y CSx^)"- 

AeA2 agA2 

Similarly, 

d 

B" = {l/2)n{n - 1)t^ Y1 E {DjkV){Djkq,)A,^ 

j,k=l mGAi 

< (i/i6)5c J2 csx^)" + E (^^^)'- 

agA2 AeA2 
By combining the above estimates we see that the left-hand side of 
(12 .61) is maj orated by 

{1-6 + Nr^) Y Qi^-^V) + J2 + I^A^n 

i/GA2 ueA AgA 

+c[2(l -5)/C( Y I^^A^I') + ((1/4)5 + iVro^) 5^ I^AV^P]. 
AeA2 AeA2 
It follows easily that by choosing tq small enough we will satisfy (12.61) 
as well as (12.51) with 5/2 in place of 6 and appropriate Ki . 

Remark 4.2. In [Hj we have seen that even Assumption 12.41 imposes 
certain nontrivial structural conditions on qx which cannot be guaran- 
teed by the size of cq if ^a is only once continuously differentiable. In 
contrast, given that Assumptions 12. l i 12.21 are satisfied and m > 2, 
we claim that Assumption 12.51 is also satisfied if cq is large enough. 
To prove our claim we notice that by (13. 5p 

AeA2 /^eA2 

so that 

n{n - 1) ^ {Sx(p)QxTx(p 
AeA2 

<m{m-l)[J2 fhvl' + NIC{ J2 M)] =■■ N'ICi M')- 

AeA2 /xeA2 /xeA2 

Now assume that c is so large that 

N' < 6c. 
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Then it follows from (12.51) that the left-hand side of fl2.6p is maj orated 
by 

(1-5) Yl Xx\Sx5M' + K^ Yl Xx\SxSM' + I, 

AgAi.i'SA^ AGAi,i/eA 

where 

/ = 2(1 - 5)clC{ J2 I^A^n + ^c/C( J2 
agA2 /^eA2 

AeA2 

We thus obtain (12.61) with 5/2 in place of 6. 

Remark 4.3. It is interesting to have sufficiently simple conditions on 
the coefficients of differential operators C which guarantee that there 
exist finite-difference schemes for which our assumptions hold. Here we 
will only give a one dimensional example. This example is based on the 
results of Remark 14.41 below, which can also be used to analyze many 
multi-dimensional situations as well in the spirit of the comments in 

Take d = 1 and 

C<^{x) = (l/2)a(x)(^"(x) + b{x)<^'{x) - c{x)<^{x). 

We assume that a > and r := ^/a, b, and c are m-times continuously 
differentiable with bounded derivatives. We take Ai = {±1} and define 

qtM = a, = {l/2)ij,b + e, 

where ^ is a constant such that > 1. By using an argument in 
Remark 6.7 of [H] and using our Remark 14. 4[ one can easily derive 
that, for a sufficiently small tq and = 1, Assumptions 12.4] and [2751 
are satisfied for all sufficiently small h (with perhaps different 6 and 
Ki) if, for n < m, 

75n'^{r'Y + 2nb' < (1 - S)c + K^a. 

Again as in [T3| we see that at points where a is close to zero either c 
should be large or b' be sufficiently negative. 

Remark 4.4. Condition (14. 2 p and its implications are discussed in many 
details in [1^ (with 2c in place of c). Here we give sufficient conditions 
for (14.21) and (14.31) to be satisfied without involving test functions ip. 
For simplicity, we only do it in case 

ta = 1 for all A G Ai. 

It is obvious that if we define = SxS^ip, then condition (14.31) can 
be rewritten in terms of What is nontrivial is that one can give 



22 



I. GYONGY AND N. KRYLOV 



sufficient conditions for fl4.3l) to hold in terms of ^\ and not the two- 
parameter object ^Xfi- In addition, we will see that these sufficient 
conditions are obtained just by slightly strengthening the corresponding 
conditions from [13] guaranteeing the ffist-order derivatives estimates. 
As in [13] one could extract further implications and simplifications 
of the new conditions of the type that on the set where c is small we 
need x\ to be uniformly bounded away from zero or px be sufficiently 
strongly monotone (see [13] for more details). 

As in Remark 14.11 we assume that Ai = — Ai and qx = q-x (> 0), 
Pa > 0, and, additionally, := ^/qx is twice continuously differentiable 
in X and is bounded on Ht along with ffist and second-order derivatives 
in X. Also we fix a constant 5 G (0, 1/4] and assume that on there 
are functions = r^x^,, Pa^ = PhXt, > 0, A, /i G Ai, such that 



m(m - l)h\6xrf,Y < 6{xx + X^^) + /^Mm' Y1 '^V < ^c, (4.5) 



By virtue of Remark 6.1 of [T3| one can always find approximations 
of the zero operator such that + L'^^ will still be approximating C 
and for the coefficients p'^ of + L'f^ we will have p'x > 1. Obviously, 
for Lh + L'f^ conditions (14.51) and (14. 6 p are satisfied with r^^ = px/j. = 
for sufficiently small h. 

For a function ^x given on Ai let us write 



and let us drop the summation sign over repeated indices in Ai. Then 
we claim that conditions (14. 2 p and (14.30 are satisfied with appropriate 
Ki, /C, and 6 if on Hj- for all functions C,x and n = 1, m we have 



AeAi 




(4.6) 



AeAi 



28^2(1 - 4:5)-\Ji + (9/2)n2(l - 4(5)~V2 + (1/2)^^3 



+25n' J2 QSxP,\ + 2nU,{SxP, + {5xr,Y) 



A,^eAi 



< (1 - 45)c|er + ^i^aXa + Sh-\x\^x + e-A 



(4.7) 



where 



/^,AeAi mSAi AgAi A,/x6Ai 
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To prove this claim, introduce 

and first recall that by Remarks 6.2 and 6.3 of after replacing there 
c with c/2 we obtain 

2n E ih^)LlTxip + 18nVi(y.) + (5/2)nV2(^) + il/2)n^Mip) 
AeAi 

< (1 - 5) E 2(^A<^) + KiQiv) + (1 - 5)c/C( E 

AeAi AgAi 

In particular condition (14.21) is satisfied. Furthermore, by substituting 
6,yip in place of ip and summing up over z/ e Ai, we get 

+ (5/2)^2 MSuV) + (1/2)^' E ^3(^.<^) 

i^GAi i/GAi 

< (1 - 5) E 2(^A^) + ^1 E Q^^'^^) + (1 - '^)^^( E i^^v^i')- 

AgA2 i^GAi AGAf 

It follows that to prove our claim, it suffices to prove that 
n{n - 1) E {SxV)QxTxV < 18^' E ^i^^-^) + (^Z^)^' E '^^(^.y') 

AGAf J^eAi ueAi 

+ (1/2)^2 J] J3(<5,y,) + (2/3)5 E 
+ (l/3)5c/C( E I'^-^^n + ^ E Q^^-*^) + E I'^-'^l')' (4.8) 

j/GAf i^GAi vGAi 

where and below by we denote generic constants independent of ip 
and {t,x) (and various e's once they appear). Observe that for A = 
(A\ A^) e Al and yU E Ai, 

Tx = l + h?5x^5x2 + /i(5ai + ^aO, = ^"'(^m + ^~^^) 
and hence 

A^Ta = + Mai(5a2((5^ + 5_^) + (5ai + ^aO(^m + '^-m)' 
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implying that 

{5\ip)Q\Txip = Si + S2, 

AGA2 

where 

S, = (1/2) Yl (^A¥^)(5Ag^)(45Ai - 

AeA2,^gAi 
AgAf ,ftGAi 

Next, as it is easy to see for A G 

SxQf, = 2{6xirf,)6x2rf, + 2rf,6xr^ 

+2h{6xirf, + 6x2r^)6xrf, + h'^i6xr^,f. 

Estimating 82- First we estimate the term S2, which contains the 
third-order differences of (f. For the main term in 5*2 we have 

Bi := 2h ^ {Sx(p){Sxirf,){6x2rf,)6xS^(p 
AeAf 

=: 16 Y M^-f) + (1/16)^, 
where by assumption (14 .Sp and Lemma 6.1 of [H] 

n{n -1)E<25J2 Qi^'^v) + ^^^/C ( ^ (5a</^)') • 

veA2 AGA2 

Hence, 

n{n - l)Bi < 16n^ ^ JiiS^^) 

+ (1/8)5 E + (l/4)(5c/C( ^(5.^)2). (4.9) 

!^eA2 AeAf 

Next, obviously, for any e > 0, (here we use that px > 0) 

B2 := 2/i Y ih^)r^,{5xrf,)5x5fj.^ 
aga2,^gAi 

AGAf./iSAi AGAf,/iGAi 

<iVs-i/C( 5^(5,(^)2) +5 5^ 2(5.^^). 

AeAi ;^ga2 
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It follows that (with e > different from the one from above but still 
arbitrary) 

n{n - 1)B2 < Ne-^K,{ i^xV?) + ^ Yl QiM- (4-10) 

Also 

AGAf ,^tGAi 

< e6c J2 i^xV? + Ne-'IC{ M') ■ 

AeA2 AgAi 

Upon combining this with (14. 9p and (14.101) we obtain 

n{n - 1)52 <{e + 5/8) ^ Q{5,ip) + Wn'' ^ M5,ip) 

i/GAJ i'SAi 

AeAf AeAi 
Estimating Si. We again start with the main term in 5*1 , which we 
split into two parts writing 

We have 

AGA2,^gAi 

= 4 ^ [ ^ iSxi5^'^))Sxr^] {S^r^)SJ^Lp 

i/,/jGAi AgAi 
u,fj,£Ai AgAi 1^,/^gAi 
i^GAi i/,/i,AGAi 

= 2 J] J2(5,(/.) + 2^ Ji(5a^). 

j^gAi AgAi 

Furthermore, 

Y ih^)ih^rf,){6x2r^)A^(p 
aga2,/^gAi 

= [ XI (^A5^(p)5Ar^] (5^r^) A^y? 

v,/iGAi AgAi 

<(V2) E [E(^A^'^^)^^^^]' + (V2) 5^ (5.r,)^(A,(p)^ 

u,ij,£Ai AgAi 1/,/iGAi 
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< (1/2) 5^ M6,v) + (1/2) J2 M^uV)- 

i^GAi j^gAi 

Next, obviously 

AGAf^jGAi AgA2 AgAi 

AgA2 AgA2 AgAi 

Finally, 

h Yl \{Sxv)Sx^S,v\ < ec6 + e-'NJC{J2i^xv)')- 

AGAf.^tGAi AgA2 AgAi 

Upon combining the above estimates we obtain 

n{n - 1)S^ < Yl [2'^i('^-^) + (5/2) J2(5.<^) + {1/2) J,{6,^)] 

i^GAi 

+ec5 Y i^x^f + ^"'^ ( E 2(5a^) + IC{Y (^a^)') ) • 

AGA2 ^AgAi AgAi ^ 

This along with (14. lip leads to (14.81) after appropriately choosing e and 
proves our claim. 
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